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ABSTRACT: We present a mathematical description of the spatial and temporal evolu- 
tion of the light intensity gradient, initiator concentration gradient, and initiation rate 
profiles for photobleaching initiator systems. This article builds on the foundation 
provided by previous researchers' contributions by generalizing the governing differ- 
ential equations to include the effects of absorption by the initiator fragments, absorp- 
tion by the monomer, and diffusion of the initiator. Simulation results have confirmed 
that, at any given time, the initiation rate profile resembles a wave front that propa- 
gates through the sample. The simulation results suggest that there is an optimum 
initiator concentration for efficient photopolymerization of thick samples. As either the 
initiator concentration or molar absorptivity is increased, the initiation rate at the peak 
of the wave front increases, the breadth of the propagating front decreases, and the rate 
of spatial propagation through the sample decreases. In contrast, the maximum pho- 
toinitiation rate and rate of spatial propagation of the initiation front can be simulta- 
neously increased with an increasing initiator quantum yield or a decreasing absorp- 
tivity of the monomer or photolysis products. Finally, on the basis of these studies, 
diffusion is expected to have negligible effects in most photobleaching systems. © 2002 
Wiley Periodicals, Inc. J Polym Sci Part. A: Polym Chem -10: 793-808, 2002; DOI 10.1002/pola. 
.10162 
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INTRODUCTION 

Photopolymerizations are generally chain poly- 
merizations in which propagating active centers 
(usually radicals or cations) are produced by a 
light-induced chemical reaction. A wide variety of 
photoinitiators are available that produce free 
radicals upon the absorption of photons of the 
appropriate energy (typically in the ultraviolet or 
visible regions of the spectrum). 1-3 The photo- 
chemical mechanism of active center generation 
has been well characterized and generally pro- 
ceeds by a cleavage, /3 cleavage, or hydrogen ab- 
straction. 4-0 In each of these reactions, the pho- 
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toinitiator is promoted to an excited state and 
then fragments (photoscission) or participates in 
an electron-transfer reaction to produce active 
centers. 

Photopolymerization has emerged as an inex- 
pensive and efficient method for producing poly- 
meric films and coatings with reduced emissions 
of volatile organic components. Light-in duced po- 
lymerization offers a number of advantages over 
the thermal curing of polymers, including spatial 
and temporal control of the initiation, high reac- 
tion rates, and lower energy requirements, which 
are important factors for a wide range of applica- 
tions. 7-9 Although radiation curing first gained 
prominence as a rapid and efficient method for 
producing films and coatings, recently, photocur- 
ing has found applications in the production of 
thick polymers and composites. 10 13 To effectively 
cure a thick section, the initiator system and 
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wavelength must be carefully selected to ensure 
that light can effectively penetrate the sample. 
Ideally, the initiation wavelength should be se- 
lected so that the initiator is the only absorbing 
species (the monomer or such should be nonab- 
sorbing). In addition, for a thick-section cure, it is 
often desirable to have a lower initiator concen- 
tration and/or extinction coefficient than for the 
photocuring of thin films. Finally, a number of 
investigators have found that it is advantageous 
to use photobieaching initiators in which light 
absorption by the initiator products is lower than 
that by the original photoinitiator molecule, 
thereby allowing more light to pass through the 
system. 

A number of initiators exhibit photobieaching 
characteristics. For example, the photolysis prod- 
ucts of benzoin ethyl ether (BEE), when irradi- 
ated with ultraviolet light in a solution of deoxy- 
genated benzene, are benzil, benzaldehyde, and 
picanol ethers.' 4 These molecules are consider- 
ably less absorbing at 365 nm than BEE itself; 
therefore, photobieaching is observed upon irra- 
diation of a solution of BEE at 365 nm. Photo- 
bleaching is particularly pronounced in solutions 
of acyl and bisacyl phosphine oxides because pho- 
tolysis results in cleavage of the carbon phospho- 
rus bond, thereby destroying the chromophore 
that absorbs in the near ultraviolet-- visible region. 
(350-410 nm). Substituted titanocenes are an- 
other class of photoinitiators that have excellent 
bleaching properties in the 440-nin region. [S 

In a photobieaching system, the light intensity 
gradient in the sample depends on both the time 
and depth beneath the illuminated surface. Ini- 
tially, the initiator concentration is uniform, and 
the light intensity will decrease exponentially 
with depth, as described by Beer's law. Immedi- 
ately after illumination, the initiator will be con- 
sumed at a rate proportional to the local light- 
intensity, thereby leading to an initiator concen- 
tration gradient with the lowest concentration on 
the exposed surface. A mathematical description 
of photoinitiation in a photobieaching system is 
complicated by the fact that the light intensity 
gradient and the initiator concentration gradient 
are coupled in this way. As recently described by 
Terrones and Pearlstein, 16 previous attempts to 
account for the variations of light intensity and 
initiation rate have failed to accurately describe 
these complex systems. Recently, two research 
groups 16,17 provided theoretical descriptions of 
initiation in photobieaching systems for the case 
in which only the original photoinitiator molecule 
absorbs (the monomer and initiator fragments 



were assumed to be nonabsorbing at the initiation 
wavelength). 

Terrones and Pearlstein 16 provided an excel- 
lent account of the shortcomings of the previous 
descriptions of photoinitiation, which did not cor- 
rectly account for the effects of the spatial varia- 
tion of the initiation rate. These authors pre- 
sented an unsteady-state, one-dimensional model 
that properly accounted for the initiator con- 
sumption and optical attenuation for free-radical 
photopolymerizations in photobieaching systems. 
Through an examination of the case in which 
only the photoinitiator absorbed at the initiat- 
ing wavelength, the resulting integro-differential 
equation could be solved analytically. These au- 
thors also effectively assumed that the sample 
was illuminated with, monochromatic light (the 
absorption of light is described with a single ab- 
sorption coefficient) and that mixing by diffusion 
or convection was negligible. Terrones and Pearl- 
stein presented results for the variations of the 
photoinitiator concentration and the photoinitia- 
tion. rate as functions of sample depth for systems 
with, an initial absorbance y (the product of the 
initiator absorption coefficient, initial initiator 
concentration, and sample thickness) ranging 
from 0.1 to 100. Interestingly, the authors found 
that with increasing initial absorbance, the local 
initiation rate becomes increasingly nonuniform 
and, for high y, assumes the form of a localized 
wave that propagates through the sample. 

In a recent contribution, Ivanov and Decker 1 ' 
also discussed the inherent inability of the stan- 
dard theoretical treatments of photopolymeriza- 
tion kinetics to account for the spatial, variation of 
the initiation rate. These authors solved the cou- 
pled differential equations that describe the spa- 
tial and temporal evolution of both the light in- 
tensity gradient and the initiator concentration 
gradient. Again, because the initiator fragments 
and the monomer were both assumed to be non- 
absorbing, the equations could be solved analyti- 
cally. The authors assumed that initiation was 
by monochromatic light, and diffusion was ne- 
glected. The initiation kinetics were coupled to 
the polymerization kinetics with a quasi-steady- 
state approximation for bimolecular radical ter- 
mination. Ivanov and Decker demonstrated that, 
under these assumptions, the light intensity gra- 
dients in the sample depended primarily on the 
absorption coefficient and initial concentration of 
the initiator. The influence of sample thickness on 
the kinetics of the photoinitiator and monomer 
consumption was investigated for systems vary- 
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ing from transparent films to optically thick sam- 
ples. 

In this article, we present a more generalized 
description of the spatial and temporal evolution 
of the light intensity gradient and the initiator 
concentration gradient in photobleaching sys- 
tems. This article builds on the foundation, pro- 
vided by Terrenes and Pearlstein 16 and Ivanov 
and Decker 1 ' by generalizing the governing dif- 
ferential equations to include the effects of ab- 
sorption by the initiator fragments, absorption by 
the monomer, and diffusion of the initiator. Be- 
cause of the complexity of the resulting set of 
coupled differential equations, a numerical solu- 
tion was obtained with the method of finite differ- 
ences. The simulation results illustrate the effects 
of a variety of system variables on the time evo- 
lution of the initiator concentration gradient, the 
light intensity gradient, and, ultimately, the pho- 
toinitiation rate profiles. 



GENERALIZED GOVERNING DIFFERENTIAL 
EQUATIONS 

For an initially homogeneous sample subject to 
uniform illumination normal to the surface, the 
evolution of the light intensity gradient through 
the depth of the sample is dependent on the con- 
centration of initiator, which also varies with 
sample depth and time. For the generalized case, 
which includes the effects of absorption by the 
initiator fragments, absorption by the monomer, 
and initiator diffusion, the following set of cou- 
pled partial differential equations describes the 
light intensity and initiator concentration pro- 
files: 

dC,(z,t) (I(z,t)\ 3 2 C-iz,t) 

dCJz,t) (I{z,t)\ d 2 CJz,t) 

¥- e 4^ c ^+ D > (2) 

- [e i C l (z,t)+A m + Sp C v (zM(z,t) (3) 

where C&,t) represents the initiator molar con- 
centration at depth z and time t; C p (z,t) repre- 
sents the photolysis product molar concentration 
at depth z and time r; I(z,t) represents the inci- 
dent light intensity at depth z and time t, with 
dimensions of energy/(area X time); e- y is the ini- 
tiator molar absorptivity with dimensions of vol- 



ume/0 ength X mole); € p is the photolysis product 
molar absorptivity, with dimensions of volume/ 
(length X mole), and accounts for the photon ab- 
sorption by all fragmentation species; N A is Avo- 
gadro's number; h is Planck's constant; v is the 
frequency of light, with dimensions of inverse 
time; &i is the quantum yield of the initiator, 
defined as the fraction of absorbed photons that 
lead to fragmentation of the initiator; D i is the 
diffusion coefficient for the initiator, with dimen- 
sions of length 2 /time; D p is the diffusion coeffi- 
cient for the photolysis products; and A m is the 
absorption coefficient of the monomer and the 
polymer repeat unit (the product of the molar 
absorptivity and the molar concentration of the 
monomer units), with dimensions of inverse 
length. 

In these governing equations, division of the 
incident light intensity I(z,t) by the product N A hv 
simply converts the intensity into the photon flux, 
with dimensions of moles of photons per unit of 
area per unit of time. Note that in eq 1, the 
change in initiator concentration with time in- 
cludes one contribution from the consumption of 
initiator in the photochemical reaction and a sec- 
ond from the diffusion of the initiator that arises 
from the resulting concentration gradient. Equa- 
tion 2, which describes the initiator products 
(fragments), has one contribution accounting for 
the production of the fragments in the photo- 
chemical reaction and a second for diffusion. Fi- 
nally, the light intensity profile through the depth 
of the sample is represented in eq 3 and includes 
terms arising from the absorbance by the initia- 
tor, monomer, and photolysis products. Note that 
any change in the monomer absorption coefficient 
upon polymerization has been neglected in the 
development of the governing differential equa- 
tions. 

The following boundary conditions and initial 
conditions apply to this system: 

CfaO) = C 0 (4) 

C.ir.Os = 0 (5) 

dC: n 

j " ------- 0 at z = 0 and 

dz 

z = z„.., x (the thickness of the sample) (6) 
/«>,*) = I 0 (7) 

Equation 4 accounts for the fact that the initi- 
ator concentration is initially uniform through 
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the depth, of the sample. Equation 5 indicates that 
photolysis product concentration is initially equal 
to zero throughout the sample, whereas eq 6 rep- 
resents the no-flux boundary condition at the 
ends of the sample for both the initiator and the 
fragments (this boundary condition is necessary 
because of the diffu.sion.al terms in the governing 
equations). Finally, eq 7 indicates that the light 
intensity at the exposed surface is equal to the 
incident light intensity at all times. 

The simplification of eqs 1-3 by the elimination 
of the diffusion term and the setting of the ab- 
sorption coefficients of the monomer arid photol- 
ysis products to zero leads to the following equa- 
tions: 



dl(z,t) 

dZ 



E^{z,i)I{z,t) 



(9) 



This simplified set of equations is equivalent to 
eqs 3 and 4 in the article by Ivanov and Decker 1 ' 
and eq 9a in the article by Terrones and Pearl- 
stein. 16 



RESULTS AND DISCUSSION 

In the remainder of this article, eqs 1-7 are solved 
simultaneously to illustrate the effects of a vari- 
ety of system variables on the time evolution of 
the initiator concentration gradient, the light in- 
tensity gradient, and, ultimately, the photoinitia- 
tion rate profiles. The simulation, conditions are 
selected to illustrate the effects of the concentra- 
tion, molar absorptivity, and quantum yield of the 
initiator, the absorptivity of the initiator frag- 
ments and the monomer, and the effect of diffu- 
sion, with the goal of identifying selection criteria 
for the effective photopolymerization of thick sys- 
tems. In future contributions, the effect of multi- 
ple wavelength initiation will be investigated. 

Representative Simulation Results 

The governing differential equations (eqs 1—7) 
were solved numerically with the method of finite 
differences to obtain profiles of the light intensity, 
initiator concentration, and initiation rates as 
functions of time and depth. For accurate results, 
a finite difference mesh size of about 1 s in the 
time domain and 0.025 cm. in the depth domain 



was used. To verify the accuracy of the numerical 
method, we compared the numerical solution to 
the analytical solution reported by Terrones and 
Pearlstein 16 for the case of no diffusion and no 
absorption by the initiator fragments and mono- 
mer. The two models gave essentially identical 
results at time zero, and any deviation would 
increase with illumination time. Therefore, the 
comparison of this model with the analytical so- 
lution of Terrones and Pearlstein 16 was made late 
in the reaction when more than 90% of the pho- 
toinitiator had been consumed. For the mesh 
sizes listed, the correlation of fit (R 2 ) for a com- 
parison of the initiator concentration profiles was 
greater than 0.994 under these stringent test con- 
ditions. This result confirms the accuracy of the 
numerical solution. 

Representative simulation results for a sample 
with a relatively low optical density are shown in 
Figures 1 and 2. These results would arise from a 
molar concentration of 0.0288 mol/L for an initi- 
ator with a molar extinction coefficient of 97.7 L 
cm" 1 mol -1 . Figure 1 contains a plot of the nor- 
malized initiator concentration [defined as C(z,i)/ 
Cy as a function of sample depth in centimeters 
for a system with no absorption by the monomer 
or the initiator fragments and no initiator diffu- 
sion. This serves as a base case to which more 
complex systems can be compared. The plot illus- 
trates that at time t = 0, the initiator concentra- 
tion is independent of depth with a uniform nor- 
malized concentration of unity. Correspondingly, 
as shown in Figure 2, the light intensity initially 
decreases exponentially in accordance with Beer's 
law. As time passes, the initiator is consumed at 
a rate proportional to the local light intensity, 
causing the concentration to decrease the fastest 
at locations near the illuminated surface. This 
results in an increased transmittance of light 
near the surface, further increasing the rate of 
consumption of the initiator. The effects of this 
bleaching process on the initiator concentration 
and the light intensity is illustrated by the family 
of curves in Figures 1 and 2. For example, note 
that the initiator concentration assumes a sigmoi- 
dal shape with nearly zero concentration near the 
illuminated surface in. Figure 1. With time, the 
region with nearly zero initiator concentration 
(and nearly 100% transmittance of light) expan ds 
as the system bleaches. For this set of simulation 
conditions, the transmittance of the sample 
reaches a normalized value of 0.9 in 84 s and a 
normalized value of 0.99 in 81 s. This relatively 
fast bleaching results from the relatively high 
light intensity, the relatively low initiator concen- 
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tration, the high initiator molar absorptivity, and 
the lack of absorption by the monomer and the 
initiator fragments. As illustrated in the next few 
sections, the initiator concentration and molar 
absorptivity both play important roles in deter- 
mining the rate of photobleaching and, therefore, 
the effectiveness of the system for thick-section, 
cure. 



Effect of the initiator Concentration 

The initiator concentration plays a key role in 
determining the time evolution of the light inten- 
sity gradient in a photobleaching s}?'stem and, 
therefore, in the effectiveness of the photopoly- 
merization of thick syst ems. A number of investi- 
gators have identified the importance of the ini- 




Figure 2. Time and depth dependence of the normalized light intensity in a 1.0- cm 
sample. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0.4 wt %), $ is 0.2, e, is 97.7 
L (mol cm) -1 , and e., is 0 L (mol cm) -1 . Mesh sizes are 1.0 s and 0.025 cm. 
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Figure 3. Effect of the various initiator concentrations on the intensity profile in a 
1.0-ern sample. The light intensity is 10 mW/mm 2 , <j> is 0.2. e f is 97.7 L (mol cm)"' 1 , and 
6p is 0 L (mol cm) -1 . Mesh sizes are 1.0 s and 0.025 cm. 



tiator concentration 10 '"' 12 ' 1 ' and have found that 
for thick-section cure it is necessary to use lower 
initiator concentrations than typically used for 
the photopolymerization of films and coatings. 
Figure 3 (a-d) illustrates the effect of the initial 
initiator concentration on the evolution of the 
light intensity gradient in a pbotobleaching sam- 
ple. This figure includes light intensity profiles for 
three different initial initiator concentrations at 
four different times, with all other simulation 
parameters being equal. Figure 3 (a) illustrates 
that, as expected, at time t = 0 (immediately 
after illumination when the initiator concentra- 
tion is still uniform), the light intensity de- 
creases exponentially with a depth of penetra- 
tion that decreases with increasing concentra- 
tion. Figure 3 (b— d) shows the light intensity 
gradients for the three initial initiator concen- 
trations at illumination times of 5, 25, and 50 
seconds, respectively. These profiles illustrate 
that the photobieaching process is more rapid 
and more efficient for the lowest initial initiator 
concentration, and that the differences in. the 



light transmittance among the three systems 
increases with increasing time. For example, 
after 50 s of illumination, the systems with 
initial initiator concentrations of 0.0268, 
0.0537, and 0.134 M exhibit a normalized inten- 
sity of 90% at depths of 0.700, 0.325, and 0.075 
cm, respectively. This effect would be even 
greater if the molar absorptivity of the initiator 
fragments were non-zero. 

Figure 4 (a-c) illustrates the effect of the ini- 
tiator concentration on the profiles of the initia- 
tion rate as a function of depth. If two active 
centers are produced upon fragmentation of the 
initiator, the instantaneous local rate of produc- 
tion of free radicals, R£z,t), is represented by the 
following equation: 16 

Rfcjt) ----- 24>s i Ciz,t)l(z,t) (10) 

With, this equation, the photoinitiation rate 
Ri(z,t) may be calculated directly from C;(z,f j and 
}{z,t. ! profiles obtained by the numerical solution 
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Figure 4. Effect of the various initiator concentrations on the photoinitiation rate in 
a 1.0-cni sample. The light intensity is 10 inW/jxini 2 , </> is 0.2. e ; is 97.7 L (mol cm)" 1 , and 
'". Mesh sizes are 1.0 s and 0.025 cm. 
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of eqs 1-7. Figure 4 (a-c) contains profiles of the 
photoinitiation rate as a function of depth for the 
same conditions shown previously in Figure 3. An 
examination of Figure 4 illust rates some interest- 
ing trends. In general, the photoinitiation rate is 
not uniform; rather, it varies sigaificasitly with 
time and depth in the sample. After the initial 
illumination, the photoinitiation rate profile ex- 
hibits a characteristic shape determined by the 
spatial variation in the product of th e light inten- 
sity and the initiator concentration. At a given 
time, with progression through the depth of the 
sample from the illuminated surface, the initia- 
tion rate increases with increasing depth to a 
maximum rate (for depths smaller than this max- 
imum, the rate is low because the initiator has 
been depleted in this region) and then decreases 
with increasing depth, because of the reduction in 
light intensity. As a result, at arry given time, the 
initiation rate profile resembles a wave front, and 
the breadth of this front is determined by factors 
such as the initial initiator concentration and the 
molar absorptivity (the breadth increases with 
decreasing initial initiator concentration and de- 



creasing molar absorptivity). This type of wave- 
like propagation through the sample was recently 
described by Terrones and Pearlstein 1 " and Ivanov 
and Decker. 1 ' 

An examination of Figure 4 (a-c) illustrates 
some trends that are important for the efficient 
photopolymerization of thick systems. For exam- 
ple, although the maximum initiation rate is 
highest for the sample with the highest initiator 
concentration, the depth of light penetration and 
the spatial propagation of the photoinitiation 
front are restricted. In contrast, although low ini- 
tiation concentrations allow efficient penetration 
of light into the sample, the photoinitiation rate 
at a given location is generally low because the 
rate is limited by the lack of initiator. These re- 
sults suggest that an intermediate optimum ini- 
tiator concentration exists for efficient photopoly- 
merization of thick samples in a timely manner. 

Effect of the Molar Absorptivity of the initiator 

The molar absorptivity of the initiator also plays 
an important role in determining the time evolu- 
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Figure 5, Effect of the initiator absorptivity on the intensity profile of a 1.0-cm 
sample. The light intensity is 10 tnW/ram 2 , C 0 is 0.0268 M (0.4 wt %), 0 is; 0.2, and e: p 
is 0 L (mol cm) -1 . Mesh sizes are 1.0 s and 0.025 cm 



tiers of the light intensity gradient in photobleach- 
ing systems. Figure 5 (a— d) contains light inten- 
sity profiles for three different initiator molar ab- 
sorptivities at four different times, with ail other 
simulation parameters being equal. Figure 5 (a) 
illustrates that, as expected, at illumination 
times near zero [Fig. 5 (a,b)], the light intensity- 
decreases exponentially with a depth of penetra- 
tion that decreases with increasing initiator ab- 
sorptivity, and the intensity at every point in- 
creases with, decreasing molar absorptivity. With 
increasing illumination time [Fig. 5 (e,d)L the sys- 
tems with higher initiator absorptivity exhibit 
higher intensity near the illuminated surface due 
to consumption of the initiator and the associated 
photobleaching. This effect arises from the fact 
that the absorban.ee of a photon is the first step in 
the bleaching process, and the system with the 
higher initiator absorptivity has a higher rate of 
photon absorption near the illuminated surface. 
Therefore, Figure 5 illustrates that the initiator 
absorptivity has two competing effects on the pho- 
toinitiation effectiveness for deep cure. A low mo- 



lar absorptivity allows more efficient penetration 
of the light into the sample; however, a higher 
molar absorptiv ity leads to higher rates of photon 
absorption and higher rates of bleaching. From 
these results, it is difficult to identify the opti- 
mum initiator absorptivity. 

For further investigation of the importance of 
the initiator absorptivity on the effectiveness of 
the photopolymerization system, it is instructive 
to examine profiles of the initiation rate as a 
function of thickness. Figure 6 illustrates the ef- 
fect of the initiator molar absorptivity on the pho- 
toinitiation rate as a function of depth for the 
same conditions shown in Figure 5. Figure 8 (a— e) 
shows that as the initiator molar absorptivity in- 
creases, the maximum initiation rate increases, 
the breadth of the propagating front decreases, 
and the rate of spatial propagation through the 
sample decreases. Figure 6 illustrates more 
clearly that the advantage of the system with the 
relatively high initiator absorptivity is a locally 
high rate of initiation that will lead to a higher 
local rate of polymerization, beginning on the sur- 
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Figure 6. Effect of the initiator absorptivity on the photoinitiation rate of a 1.0-cm 
sample. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0,4 wt %), 4> is 0.2, and e p 
is 0 L (mo! cm) -1 . Mesh sizes are 1.0 s and. 0.025 cm. 



face. This could allow a barrier polymer to form on. 
the surface, whereas the relatively sharp initia- 
tion front mo ves through the depth of the sample. 
The advantage of the system with the low absorp- 
tivity results from the more efficient light, pene- 
tration into the sample. These systems exhibit 
appreciable photoinitiation rates deep into the 
sample in a relatively short period of time. This 
result is important for systems that require the 
development of polymer properties soon after il- 
lumination. An initiation system exhibiting both 
of these advantages could be achieved with two 
different incident initiation, wavelengths with one 
or more initiators. 

Figures 3-6 illustrate the importance of exam- 
ining the effects of the initiator concentration and 
absorptivity separately. Because the product of 
the initial initiator concentration and the initia- 
tor absorptivity is proportional to the initial opti- 
cal density, it is tempting to characterize these 
two parameters with a single parameter. How- 
ever, it is important to optimize these two param- 
eters separately because each has distinct effects 



on the time evolution of the light intensity gradi- 
ent and the initiation rate profile. For example, 
consider two systems of equal thickness with the 
same optical density: one with a high initial ini- 
tiator concentration and a. low absorptivity and 
the second with a low initial initiator concentra- 
tion and a high absorptivity. Immediately upon 
illumination, the two systems will exhibit identi- 
cal light intensity profiles throughout the sample 
because the initial optical densities are the same. 
With increasing time, however, the initiator con- 
centration will change because of photobleaching, 
whereas the absorptivity remains constant. The 
system with the higher initiator absorptivity will 
photobleach more rapidly because a higher frac- 
tion of the molecules absorb photons, and the 
photoinitiation front will propagate more rapidly 
through the thickness of the sample. 

Effect of the initiator Quantum Yield 

It is interesting to contrast the effect of the initiator 
quantum yield with the effect of the initiator molar 
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Figure 7. Effect of the initiator quantum yield on the photoinitiation rate of a 1.0- cm 
sample. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0,4 wt %), e t is 97.7 L (mol 
cm) -1 , and e is 0 L (mol cm)" 1 . Mesh sizes are 1.0 s and 0.025 cm. 



absorptivity discussed previously. Although the mo- 
lar absorptivity characterizes the probability that 
the initiator will absorb an incident photon, the 
quantum yield describes what happens after the 
photon is absorbed. We assume that the initiator 
molecules that are not dissociated after absorption 
will return to the ground state by either radiative or 
nonradiative processes and will be available for 
further photochemical reactions (the recent arti- 
cles by Terrenes and Pearlstein 56 and by Ivanov 
and Decker 17 are based on the same assumption). 
Figure 7 (a,b) illustrates the effect of the initiator 
quantum yield on the photoinitiation rate as a 
function of sample depth. The figure illustrates 
that as the quantum yield increases from 0.2 to 
0.4 and finally to 0.8 with all other simulation 
variables held constant, the maximum photoini- 
tiation rate increases and the rate of spatial prop- 
agation of the polymerization front increases. 
This is in contrast to the effects of both the initial 
initiator concentration and the initiator absorp- 
tivity. For both of those variables, an increase in 
the maximum photoinitiation rate had an associ- 
ated decrease in the rate of propagation of the 
photoinitiation front. 

Effect of the Photofysss Product Absorptivity 

The effect of the photolysis product absorptivity 
on the time evolution of the light intensity gradi- 
ents is illustrated in Figure 8 (a-d). This figure 
contains light intensity profiles at four different 
times for three different molar absorptivities of 
the photolysis products: e p - 0, which corresponds 
to perfect bleaching; e p = 6.5 L/(mol cm), which 
corresponds to an. intermediate level of bleaching, 



and e p = 97.7 L/(mol cm), which is the same 
absorptivity as the original initiator and, there- 
fore, corresponds to a completely nonbleaching 
system. Immediately upon illumination, the three 
systems exhibit identical light intensity profiles 
because there are no photolysis products at time 
zero. With increasing time, the light intensity 
gradient for the system without bleaching re- 
mains unchanged, whereas the photobleaching 
systems allow more efficient light penetration 
into the sample. Figure 9 (a— c) illustrates the 
effect of the photolysis product molar absorptivity 
on the photoinitiation rate profile for the same 
conditions shown in Figure 8. Figure 9 illustrates 
that for the nonbleaching system, because the 
light intensity gradient is stagnant, the photoini- 
tiation rate remains close to zero at the unillumi- 
nated end of the sample, and the photoinitiation 
rate at a given point decreases with time as the 
photoinitiator is consumed. Therefore, it is diffi- 
cult to achieve thick-section cure for a nonbleach- 
ing photoinitiator, and this effect is magnified for 
systems as the initiator concentration or sample 
thickness is increased. Figure 9 illustrates that as 
the molar absorptivity of the photolysis products 
is reduced, the rate of propagation of the photo- 
initiation front increases, and for the perfectly 
bleaching case, the peak in the photoinitiation 
rate reaches the bottom of the sample within 50 s 
for the conditions in the simulation. This time 
could be reduced further with decreasing initiator 
concentration or increasing incident light inten- 
sity. Again, these results illustrate the impor- 
tance of proper system design for the efficient 
photo-polymerization of thick systems. 
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Figure 8. Effect of the photolysis product absorptivity on the intensity profile of a 
1.0-cm sample. The Light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0.4 wt %), <j> is 0.2, 
and 6; is 97.7 L (moi cm) " 1 . Mesh sizes are 1.0 s and 0.025 cm. 



Effect of the Monomer Absorption Coefficient 

The effect of the monomer absorption coefficient 
on the time evolution of the light intensity gradi- 
ents is illustrated in Figure 10 (a— d). This figure 
contains light intensity profiles at four different 
times for three different absorption coefficients of 
the monomer/polymer. For convenience, the ab- 
sorption coefficient (which is a product of the mo- 
lar absorptivity and the molar concentration) was 
chosen to characterize the tendency of the mono- 
mer/polymer to absorb photons. Even if the molar 
absorptivity of the monomer is very low, the ab- 
sorption coefficient of the monomer may be signif- 
icant because the monomer molar concentration 
is generally orders of magnitude higher than that 
of the initiator. Again, any change in the absorp- 
tion coefficient upon polymerization has been ne- 
glected in the development of the governing dif- 
ferential equations. 

For the construction, of Figure 10, three repre- 
sentative values of the monomer absorption coef- 
ficient have been selected: A m ------- 0 cm™ 1 , which 



corresponds to a transparent matrix; A m = 0.3 
cm" 1 ; and A m - 6 cm -1 . The latter value corre- 
sponds to the case in which the monomer absorp- 
tion coefficient matches that of the initiator at 
time zero. This figure clearly illustrates the im- 
portance of selecting an initiation wavelength for 
which the monomer is transparent. Although the 
two systems with relatively low absorption by the 
monomer allow the light to penetrate the sample 
with increasing illumination time, the system 
with the high monomer absorption coefficient has 
a relatively stagnant light intensity gradient. The 
effect of monomer absorption on the photoinitia- 
tion rate profile is shown in Figure 11. The system 
with high absorption of light by the monomer 
exhibits both a lower maximum initiation rate 
and a relatively stagnant photoinitiation front 
that reaches only 20% of the sample in the time it 
takes for the other systems to propagate com- 
pletely through the depth of the sample. Figure 
11 makes it clear that the absorption of light by 
the monomer/polymer makes much less efficient 
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Figure 9. Effect of the photolysis product absorptivity on the photoinitiation rate of 
a I.O-cm sample. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0.4 wt %), 6 is 0.2, 
and 6j is 97.7 L (moi cm) " 1 . Mesh sizes are 1.0 s and 0.025 cm. 



use of the incident photons because the integral of 
the photoinitiation rate profile is greatly reduced. 
Therefore, by selecting an initiation wavelength 
that the monomer does not absorb, we can signif- 
icantly increase the maximum photoinitiation 
rate (indeed, it is possible to increase the photo- 
initiation rate throughout the sample) and in- 
crease the rate of spatial propagation of the poly- 
merization, front. All of these factors lead to more 
efficient photopolymerization of thick systems. 

Effect of the Initiator Diffusion 

Diffusion of the unreacted initiator and the pho- 
tolysis products could arise from the concentra- 
tion gradients produced in the photobleaching 
process. The relative importance of diffusion will 
depend on the value of the diffusion coefficient, 
the concentration gradient, and the time required 
for photobleaching to be complete. Diffusion will 
be more important for photobleaching systems 



containing high initiator concentrations, which 
lead to steep concentration gradients and rela- 
tively long times for photobleaching. To explore 
the importance of diffusion in a typical thick- 
section cure, we have examined the two limiting 
cases of no diffusion and perfect mixing (which 
corresponds to an infinitely large diffusion coeffi- 
cient). Figure 12 contains a plot of the normalized 
initiator concentration as a function of depth for a 
representative photobleaching system after 25 s 
of illumination. The figure contains three curves: 
one for perfect mixing (which is a horizontal line 
because the concentration is uniform in this case), 
one for a representative diffusion coefficient of 
0.001 cm 2 /s (we chose a relatively large value of 
the diffusion coefficient to exaggerate the diffu- 
sional effects), and a third corresponding to the 
lack of diffusion. Similarly. Figure 13 contains a 
plot of the light intensity as a function of the 
sample depth for the same system shown in Fig- 
ure 12. The system corresponding to perfect mix- 
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Figure 10. Effect of the monomer absorption coefficient on the intensity profile of a 
1.0-ern sample. The light intensity is 10 niW/min 2 , C 0 is 0.0268 M (0.4 wt %), <j> is 0.2, 
gj is 97.7 L (mol cm)" 1 , and e p is 0 L (mol cm) -1 . Mesh sizes are 1.0 s and 0.025 cm. 



ing exhibits an exponential decay in accordance 
with Beer's law, whereas the other two systems 
exhibit the sigmoidal shape characteristic of pho- 
tobleaching systems. These results suggest that 
diffusion plays at most a moderate role in these 
photobleaching systems. In both Figures 12 and 
13, the representative system with substantial 
diffusion, lies much closer to the limiting system 
with no diffusion than to the limiting case of perfect- 
mixing. Furthermore, in most photopolymerization 
systems, the operative diffusion coefficient is likely 
to be at least an order of magnitude smaller than 
the value of 0.001 cm% that we chose for these 
studies. In addition, the diffusional mobility of the 
photoinitiator will be further restricted by the for- 
mation of the polymer matrix in these photopoly- 
merization systems. On the basis of these studies, 
we expect diffusion to have negligible effects in most 
photobleaching systems because the cliffusional 
timescale is large relative to that of pbotoinitiation. 
Diffusion will play the largest role in systems con- 
taining high photoinitiator concentrations and will 
tend to broaden the propagating photoimtiatioii 
front in these systems. 



CONCLUSIONS 

In this article, we have presented a description of 
the spatial and temporal evolution of the light 
intensity gradient, initiator concentration gradi- 
ent, and initiation rate profiles for photobleaching 
initiator systems. This article builds on the foun- 
dation provided by researchers' previous contri- 
butions by generalizing the governing differential 
equations to include the effects of absorption by 
the initiator fragments, absorption by the mono- 
mer, and diffusion of the initiator. Because of the 
complexity of the resulting set of coupled differ- 
ential equations, a numerical solution has been 
obtained with the method of finite differences. 
Simulation results have confirmed that, at any 
given time, the initiation rate profile resembles a 
wave front, and the breadth of this front is deter- 
mined by factors such as the initial initiator con- 
centration and the molar absorptivity. 

The simulation results suggest that there is an 
optimum initiator concentration for the efficient 
photopolymerization of thick samples. For initia- 
tor concentrations above the optimum, the depth 
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Figure 11. Effect of the monomer absorption coefficient on the photoinitiation rate of 
a 1.0-cm sample. The light intensity is 10 mWAnm 2 , C 0 is 0.0268 M (0.4 wt %), tj> is 0.2, 
q is 97.7 L (mol cm)" and g is 0 L (mol cm)'" 1 . Mesh sizes are 1.0 s and 0.025 cm. 



of light penetration and the spatial propagation of 
the photoinitiation front are restricted, whereas 
for concentrations below the optimum, the photo- 
initiation rate is limited by the lack of initiator. 
As the initiator concentration is increased, the 
initiation rate at the peak of the wave front is 
increased, although the rate of propagation of th e 
front through the sample is decreased. 

The initiator absorptivity has two competing 
effects on the photoinitiation effectiveness for 
deep cure. A low molar absorptivity allows more 
efficient penetration of light into the sample; how- 
ever, a higher molar absorptivity leads to higher 
rates of photon absorption and higher rates of 
bleaching. In addition, as the initiator molar ab- 
sorptivity increases, the maximum initiation rate 
increases, the breadth, of the propagating front 
decreases, and the rate of spatial propagation 
through the sample decreases. 

Simulation results have revealed that it is im- 
portant to optimize the initiator concentration 
and absorptivity separately. For example, for two 
systems with the same optical density, the system 



with the higher initiator absorptivity will photo- 
bleach more rapidly, and the photoinitiation front 
will propagate more rapidly through, the thick- 
ness of the sample. In contrast to the effects of the 
initial initiator concentration, and absorptivity, 
an increase in the initiator quantum yield results 
in an increase in both the maximum photoinitia- 
tion rate and the rate of spatial propag ation of the 
polymerization front. 

The simulation results illustrate the impor- 
tance of photobieaching for effective thick-section 
cure. As the molar absorptivity of the photolysis 
products is reduced, the rate of propagation of the 
photoinitiation front increases (the light intensity 
gradient and, therefore, the photoinitiation front 
are completely stagnant in the nonbleaching 
case). Under simulation conditions, the photoini- 
tiation front can reach a depth of 1 cm in less than 
a minute. The simulation, results have confirmed 
that the absorption of light by the monomer or 
polymer makes much less efficient use of the in- 
cident photons. Therefore, by selecting an initia- 
tion wavelength that the monomer does not ab- 
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Figure 12. Effect of the initiator diffusion on the concentration profile of a 1.0-cm 
sample at 25 s. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0.4 wt %), <j> is 0.2, 
6j is 97.7 L (mol cm)" 1 , and « is 0 L (mol cm)" 1 . Mesh sizes are 1.0 s and 0.025 cm. 



sorb, we can significantly increase the maxi- 
mum photo initiation rate (indeed, it is possible 
to increase the photoinitiation rate throughout 
the sample) and increase the rate of spatial 
propagation of the polymerization front. Fi- 
nally, on the basis of these studies, we expect 
diffusion to have negligible effects in most pho- 



tobleaching systems because the diffusional 
timescale is large relative to that of photoinitia- 
tion. Diffusion will play the largest role in sys- 
tems containing high photoinitiator concentra- 
tions and will tend to broaden the propagating 
photoinitiation front in these systems. This 
model provides a valuable tool for evaluating 




Figure 13. Effect of the initiator diffusion on the intensi ty profile of a 1.0-cm sample 
at 25 s. The light intensity is 10 mW/mm 2 , C 0 is 0.0268 M (0.4 wt %), 6 is 0.2, e s is 97.7 
L (mol cm) -1 , and e., is 0 L (mol cm) -1 . Mesh sizes are 1.0 s and 0.025 cm. 
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the effectiveness of photobleaching initiator 
systems for thick-section cure. 

This article is based on work supported by the Indus- 
try/University Cooperative Research Center for Funda- 
mentals and Applications of Photopoiymerization. 
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